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Abstract 


DIFFUSION  OF  A  PLASMA  ACROSS  A  MAGNETIC  FIELD 

by 

Ronald  Stephen  Finn 


An  investigation  of  the  problem  of  plasma  diffusion  through 
a  magnetic  field  is  made  through  the  macroscopic  fluid  equations. 
The  problem  is  set  up  in  a  simple  one-dimensional  geometry  of  a 
semi-infinite  half-space  filled  with  plasma  facing  a  semi-infinite 
half-space  containing  a  uniform  magnetic  field,  and  the  equations 
are  used  to  describe  their  subsequent  behavior.  An  attempt  is 
made  to  solve  the  equations  analytically,  but  their  final  solution 
is  obtained  numerically.  The  solutions  are  also  used  to  determine 
temperature  and  entropy  increases  in  the  plasma. 


DIFFUSION  OF  A  PLASMA  ACROSS  A  MAGNETIC  FIELD 


by 

Ronald  Stephen  Finn 


I  Introduction 

There  has  been  very  recently  a  greatly  Increased  Interest  in 
scientific  and  engineering  applications  of  ionized  gases  and  plasmas; 
this  has  demanded  from  known  theory  a  quantitative  idea  of  how  plasmas 
behave  under  given  conditions.  An  example  of  these  is  the  problem  of 
plasma  containment  by  a  magnetic  field;  or^  stating  it  more  generally, 
the  interaction  of  a  plasma  with  a  magnetic  field  generated  by  ex¬ 
ternal  sources . 

This  plasma-field  interaction,  which  is  the  subject  of  this  dis¬ 
cussion,  has  a  number  of  aspects  that  complicate  its  analysis;  the 
complications  resulting  in  non-linear  dlfferentlEd.  equations.  As  a 
consequence,  the  investigation  of  this  interaction  is  limited  to  a 
simple  example;  this  is  to  make  the  equations  concerned  tractable 
enough  to  yield  quantitative  data  with  a  reasonable  amount  of  effort. 
To  illustrate,  the  processes  considered  can  be  crudely  described  as 
"slow"  motions  of  plasmas  through  "strong"  magnetic  fields,  in  which 
the  electromagnetic  and  thermodynamic  properties  of  the  plasmas  do 
not  change  appreciably.  What  is  meant  by  "slow"  motions  and  "strong" 
fields  will  be  made  clearer  further  on.  In  addition  to  this,  the 
geometry  of  the  problem  is  reduced  to  its  simplest  form.  Doing  this 
naturally  limits  the  applicability  of  the  resulting  data;  still,  it 
produces  a  good  approximation  of  the  plasma-field  Interaction  in  a 
variety  of  cases. 


II  General  Discussion 
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The  formal  statement  of  this  problem  most  amenable  to  quali¬ 
tative  investigation  is  in  terms  of  the  so-called  macroscopic 
equations  of  motion.  Essentially  these  are  the  equations  that 
result  from  averaging  the  motions  of  many  particles;  this  averaging 
being  done  by  integration  of  the  Boltzmann  equation  in  velocity 
space.  By  means  of  this  integration  Spitzer^^)  defines  a  current 
density  j,  an  average  fluid  velocity  u,  and  derives  an  equation  of 
fluid  motion  and  a  genercdized  form  of  Ohm's  law. 

In  the  case  to  be  considered  here,  the  fluid  is  assumed  to  be 
composed  only  of  singly  charged  ions  of  identical  masses  and  of 
electrons.  To  generalize  to  a  number  of  types  of  ions  does  not 
complicate  the  situation,  but  it  turns  out  that  this  does  not  lead 
to  any  additional  information.  A  monatomic  gas  was  chosen  simply 
to  make  the  ratio  of  specific  heats  equal  to  ^/3  in  calculations 
done  later.  For  this  case,  then,  the  equation  of  motion  takes  the 
form: 

+  (0(u  .  7)u  -  j  X  B  -  7p  (1) 

where  B  is  the  magnetic  field  strength,  p  is  the  scalar  pressure, 
euid  (O  is  the  mass  density  of  the  ions.  This  is  simply  Euler's 
equation  with  the  addition  of  a  j  x  B  force. 

There  are  certain  assumptions  already  inq)liclt  in  this  equation: 
the  mass  of  the  electrons  and  the  gravitationeO.  force  are  both  neg¬ 
lected,  and  the  number  densities  of  the  electrons  and  ions  are  taken 
to  be  exactly  equal  at  all  points  in  the  fluid.  Also,  viscosity 
effects  are  neglected,  allowing  what  is  actuedly  the  stress  tensor 
to  be  taken  as  a  scalar  pressure.  Accounting  for  directional  pres¬ 
sure  differences  due  to  external  magnetic  fields  by  using  a  stress 
tensor  with  at  least  different  non-zero  diagonal  terras  might  seen 
actually  necessary  for  this  type  of  problem;  however,  it  will  turn 
out  that  this  does  not  make  any  difference  in  the  calculations. 

Also  for  this  case,  the  generalized  Ohm's  law  reduces  to  the 

form: 

■  E  +  u  X  B  (2) 


where  ^  is  the  electrical  resistivity  and  E  is  the  electric  field. 
One  assumption  made  here  is  that  the  resistivity  is  a  constant  scalar; 
once  again,  writing  it  as  a  diagonal  tensor  to  account  for  external 
magnetic  fields  does  not  affect  the  calculations.  It  must  be  admitted 
that  neglecting  its  dependence  on  the  temperature  considerably  limits 
the  usefulness  of  the  solution;  nevertheless,  this  is  the  first  as¬ 
sumption  that  will  have  to  be  made  strictly  for  simplification. 

Another  assunqption  is  that  a  ^  j/  ^  t  term  is  negligible  small  com¬ 
pared  to  the  others;  here  the  limitation  of  "slowness"  is  first 
applied. 

The  two  equations  above  resulted  from  the  "plasma"  properties 
of  the  fluid  in  this  problem.  The  remaining  equations  to  be  intro¬ 
duced  are  of  a  more  general  nature.  Firstly,  there  are  the  continuity 
equation 

+  V  .  (  ^  5)  -  0  (3) 

and  the  equation  of  state 

p  -  ^  R  T  (4) 

where  R  is  the  gas  constant  and  T  is  the  teDq>erature .  In  addition, 
there  are  the  Maxwell  equations 

V  X  E  -  -  -^  (5) 

and 

j  --i-  V  X  B  (6) 

The  currents  due  to  polarization  and  the  time  rate  of  change  of  the 
electric  field  have  been  neglected  as  being  small  conqpared  to  the 
magnetization  current  which  is  Included  in  the  V  x  B  term  in  this 
last  equation.  The  permeability,  is  taken  equal  to  that  of  free 
space. 

Finally,  there  is  the  energy  equation 
cv(-^  +  u  .  7)  ln(pp"’'  )  »  ^4" 

where  Cy  is  the  constant-volume  heat  capacity  and  Y*  is  the  ratio  of 
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specific  heats.  This  equation  iti5>lie8  two  very  important  assumptions j 
one,  that  the  only  source  of  energy  is  ohmic  heating  (viscosity  having 
already  been  neglected);  and  the  other,  that  there  is  no  heat  transfer 
in  the  fluid,  i.e.,  all  processes  are  adiabatic.  This  last  assumption 
Is  valid  when  heat  diffusion  is  ouch  slower  than  the  process  to  be 
described. 

Now  tabulating  the  equations  obtained: 


f  41  *  P  ■  j 

(1) 

'•JJ-E+uxB 

(2) 

+  K  (p  S)  -0 

(3) 

p  -  p  R  T 

(h) 

(5) 

3  ?  X  B 

>*0 

(6) 

Cv(-^  +  u  .  7)  ln(pp  )  « 

(7) 

it  is  seen  that  there  are  seven  equations  containing  seven  variables. 

It  may  also  be  seen  that  the  equations  are  non-linear,  which  imnedlately 
precludes  any  straightforward  analytical  Integration  of  them.  Thus, 
proceeding  with  their  investigation  requires  some  additional  simpli¬ 
fications  . 

% 

Before  continuing  with  this,  it  Is  instructive  to  roughly  sketch 
some  properties  of  plasmas  in  strong  magnetic  fields  produced  by  con¬ 
stant  external  currents.  Here  "strong  magnetic  fields"  means  fields 
that  are  much  larger  than  any  field  set  up  by  currents  in  the  plasma. 

To  the  first  order,  the  magnetic  force  lines  tend  to  be  "frozen"  in 
a  highly  conducting  gas.  If  a  motion  is  inparted  to  the  gas,  the 
force  lines  are  "dragged  along"  with  It.  To  the  next  order,  current 
dissipation  due  to  the  finite  resistivity  of  the  gas  allows  a  relative 
motion  between  the  gas  and  the  field  lines.  To  see  this  more  clearly, 
one  can  imagine  a  cluster  of  externally  produced  parallel  magnetic 
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field  lines  in  a  plasma  with  a  finite  resistivity.  Because  the 
plasma  is  resistive,  the  current  around  an  area  in  it  will  have  to 
decay  if  it  is  not  regenerated  by  externally  intensifying  the  fields. 
Thus  the  number  of  flux  lines  passing  through  the  area  will  have  to 
be  reduced,  and  if  their  total  number  is  constant,  some  of  them  must 
appear  outside  of  the  area.  Therefore,  the  cluster  spreads  out  and 
trtiat  the  whole  process  "looks  like"  is  a  diffusion  of  the  field 
lines  through  the  plasma. 


Ill  The  Infinite  One-PimensioneLl  Case 

Keeping  in  mind  the  limitations  of  "strong"  magnetic  fields  and 
"slow"  processes,  it  is  now  possible  to  set  up  a  specific  example  of 
this  problem  that  will  lend  Itself  to  further  simplification  of  the 
seven  equations.  Consider  at  time  t  ■  0  a  perfectly  uniform,  com¬ 
pletely  ionized,  monatomic  plasma  of  density  and  resistivity 
extending  from  the  yz-plane  to  infinity  along  the  positive  x-axls. 

The  plasma  is  not  allowed  to  have  any  motion  which  is  not  along  the 
x-axis;  i.e.,  its  motion  is  one-dimensional.  Consider  also  a  per¬ 
fectly  \iniform  magnetic  field  of  intensity  Bq  at  t  o  o  extending  from 
the  yz-plane  to  infinity  along  the  negative  x-axls.  The  field  lines 
are  all  perfectly  straight  and  parallel  to  the  z-axls,  pointed  in  the 
positive  z-dlrectlon,  and  extend  to  infinity  in  both  z-dlrectlons 
(Figure  l). 


Figure  1 


Because  the  field  Is  infinite  but  does  not  extend  past  x  ■  0,  there 
is  at  t  0  an  infinite  current  sheet  with  the  current  travelling  in 
the  negative  y-direction  and  covering  the  x  ■  0  surface  of  the  plasma. 

To  seek  a  solution  for  the  behavior  of  this  system  for  times 
t  >  0,  one  first  considers  the  equation  of  motion 

P  +  P(u  .  7)u  -  J  X  B  -  7p  (l) 

Since  the  plasma  is  always  uniform  in  any  plane  perpendicular  to  the 
x-axls^  the  V-operator  reduces  to  Vi  X.  Also,  if  the  B-field  is 
always  considered  in  the  z-dlrection,  the  net  currents  in  the  plasma 
always  flow  in  the  y-direction.  Then  the  equation  of  motion  becomes 


How  it  may  be  seen  why  it  was  not  necessary  to  regard  the  resistivity 
and  the  pressure  as  tensors.  Even  if  the  tensor  form  had  been  re¬ 
tained,  only  one  term  in  each  would  have  been  used. 

At  this  point  the  criterion  of  "slowness"  is  again  applied  to  the 
motion  of  the  plasma.  This  means  tha,t  the  terms  on  the  left  are  neg¬ 
ligibly  small,  or  that  the  terms  on  the  right  nearly  cancel.  Thus, 

JB  .  p  (9) 

Since  for  this  configuration  the  equation 

J  .  -i-  f  X  B  (6) 

>^0 


becomes 


.3  -  - 


1  ^B 


(10) 


because  B  only  has  a  functional  dex>endence  on  the  x-coordinate,  the 
equation  of  motion  then  becomes 


(11) 


nnri  this  may  be  integrated  to  give; 
b2 


P  + 


constant 


(12) 


Thus,  b2/^^  may  be  considered  as  a  "magnetic  pressure",  and  the 
total  pressure  at  aU  points  in  space  remains  constant  as  the  magnetic 
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and  plasma  pressures  vary.  Simply  from  this  equation  it  is  possible 
to  see  that  the  plasma  leaks  slowly  into  the  magnetic  field  while  the 
magnetic  field  correspondingly  leaks  into  the  plasma.  Neither  the 
plasma  nor  the  field  has  a  sharp  boundary  at  t  >  0;  and  the  current 
setting  up  the  magnetic  field  becomes  distributed  through  the  plasma 
rather  than  being  concentrated  in  a  sheet.  A  plot  of  plasma  density 
and  field  strength  might  then  be  expected  to  look  like  Figure  2. 


Figure  2. 


The  next  series  of  steps  in  seeking  a  solution  to  the  equations  consists 
of  reducing  them  as  far  as  possible.  The  calcxilations  are  first  simpli¬ 
fied  by  using  the  facts  that  the  magnetic  fields  are  only  the  the  z~ 
direction  and  the  currents  and  therefore  electric  fields  are  only  in 
the  y-direction  (there  are  no  external  fields);  and  that  all  quantities 
depend  functionally  on  the  x-dlrection.  Ohm's  lav  then  becomes 


J  «  E  -  u  B 

and  differentiating  it  with  respect  to  x  yields 


The  continuity  equation  becomes 


and  Faraday's  law  becomes 


(13) 

(14) 

(15) 


ix  “  "  it 


(16) 


-  8  - 


Substltutlxig  for  ^u/^x  from  the  continuity  equation,  and  for  ^E/^x 
from  Faraday's  law  into  the  differentiated  form  of  Ohm's  law  yields 


4l.„is*l(44,„4£) 


d  t  “  ^  X  ^  ^  d  X 

Then  substituting  for  j  from  the  first  Maxwell  equation  yields 


(17) 


r  ■  <-5t  *  “  -5^)  -  jt  *  “  (18) 


Now  the  process  of  reduction  is  carried  through  starting  from  the 
energy  equation 

p  T 

Substituting  for  p  T  from  the  equation  of  state  and  dividing  by  Cy 
yields 


"  -It) 


(7) 


1/dp  .  dp\  Y’/dP  .  dPs 

;(-n  *  “  T?'  ■  f  <TE  *  “ 


(19) 


1  ^  B 

Multiplying  by  p  and  substituting  R/d^-  Y^-  1  and  J  *  -  ^ 


Po' 

Now  the  e3q)resBion  for  p  in  terms  of  B  derived  from  the  equation  of 
OK>tlon  is  substituted,  yielding 

B  /SB  dB^  _  V  -  1/ dp  .  .  dPN  *1  ,^/dB^2 

-  ji^Tt  * « -ji;)  -  r-^i!(Tt  *  “  ix-r;) 

(21) 

Simple  substitution  produces  no  further  reduction  after  this  point, 
and  the  process  must  be  continued  by  the  following  transformation  of 
coordinates :  let 


dx 


T  -  t 


(22) 


Then 


and 


d  dx  d  dT  d 

^x"  d*  d^  dx  dT 

JLA- 

"  Po 

^  dx  d  ,  dT  d 

t  *  dt  dX  Tt  dT 


(23) 

(24) 
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l.iut 


“it  /"St®*  *  ■  "t 

/•«  ^rnm 


p. 


by  the  equation. of  continuity,  so  that 

^  Pu  ^  ^ 

Tt-'-^-Sx^Tr 


Thus 

— i^"5x^iT^''p“"5x’i;T 

and  substituting  appropriately  in  Eqs.  (l8)  and  (2l)  yields 


“*1  d  /  ifx 

M;  8  TX^8  "SX^ 


if  fig 

It  ■  g  Tr 


-f 


^f  ^i-f^ 

^T  ■  *  2 


i  Ae 

g  iT  M* 

using  the  dimensionless  variables  g 


(r 

p  /p. 


iX' 
;  f 


(25) 


(26) 


(27) 


(28) 


(29) 


B/Bc 


This  transformation  has  served  its  intent  by  eliminating  u  as  a 
variable;  but  it  has  also  had  other  effects  in  changing  the  x-axls. 

The  new  axis  is  "squeezed"  compared  to  the  real  x-axis  when  the  pro¬ 
portionality  factor  P  / P*  is  small.  Also,  X  is  never  less  than  zero 
because  the  defining  integral  is  never  negative.  Thus,  large  negative 
values  on  the  x-axis  correspond  to  small  positive  values  on  the  X-axis, 
while  large  positive  values  on  the  x-axis  correspond  to  large  positive 
values  on  the  X-axis.  However,  this  transformation  has  reduced  the 
problem  to  two  equations  in  two  lonknowns,  and  the  remaining  discussion 
is  involved  with  trying  to  solve  them. 

•For  further  convenience,  the  two  equations  will  be  rewritten  in 
terms  of  one  variable 


,  E  JL 
V  n  vs 


1  vt 

not  to  be  confused  with  the  z-coordinate,  so  that 


■5x 


■5xTI  “Vt  ;FTi  -xT^ 


and 


T  Vt 

i  /jS  JL  -i_  ■  -1 1 

fr  “  iT  T  ^  ®  T  i  z 


(30) 


(31) 


(32) 
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Eq.  (28)  then  becomes 

z  ^  /  z  if  X  _  1  z  i  f  .  1  f  Z  ig 

M*  ®  X  i  z'®  X  ‘^T  dz  ^gT 

and  expanding  the  term  on  the  left  yields 

"(  z  r  z  i  g  i  f  z  i^f  g  i  f  z  i  f  ix 

-:i8xLxizT^^x8  7;2-"xTI'®;^TST^ 

+  iii  +  fi4^ .  fi  4-^ 
>*.  ®x[xirzTz  X®^22^xTz  xTz 

A  ^  fl£  ll  +  a 
=  -  i».S2[^^  ®  iz^J 

Thus  Eq.  (33)  becomes 

z  I^Tz  ^  ^  J  s  i  Z 

Correspondingly,  Eq.  (29)  becomes 

2  T  ^  ^  T  g  T  7z  ^  ^  i  z' 

or  5 

Thus,  writing  primes  for  i/dz-s 
•^(g'f  +  gf")  -  -f  + 


-  g' 


ff  +|(i-f2)ig'  -2(r- 1)  f-(f)‘ 


(33) 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


It  may  Immediately  be  shown  that  these  equations  lead  to  a  pro¬ 
cess  that  at  least  resembles  ordinary  diffusion.  Taking  the  limiting 
case  of  an  incooqpresslble  plasma;  i.e.,  g  — ►  1,  g*  ^  0,  and  Y  ►  00 

causes  Eq.  (39)  to  yield  the  value  fol  and  Eq.  (38)  to  become 


2 


-f 


Keurranging  this  as 

=  (In  f')'  - 
and  integrating  it  yields 

-  -  • 

"  F  ' 


z 

2 


(1^0) 

(41) 


■z^/k 


In  f 


(42) 
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ejid  integrating  it  again  yields 

s 

e  *  dz  (43) 

Thus  f  in  this  case  is  of  the  form  of  the  error  function.  What  this 
means  is  that  the  plasma  sits  perfectly  still,  being  incompressible 
and  infinite,  and  the  magnetic  field  diffuses  into  it.  This  is  exactly 
the  process  undergone  by  a  magnetic  field  diffusing  into  a  motionless 
meteil  conductor. 


IV  Analytical  Integration 


Because  Eqs.  (38)  and  (39)  are  non-linear  second  order  differential 
equations,  they  cannot  be  straightforwardly  integrated  over  all  values 
of  z.  What  can  be  done,  however,  is  to  look  at  the  forms  of  these  equa¬ 
tions  as  f  and  g  asymptotically  approach  values  determined  from  physical 
reasoning.  Remembering  that  the  functions  are  now  in  terms  of  the  trans¬ 
formed  coordinate  X,  it  may  be  seen  that  at  large  z,  g  approaches  one 
and  f  approaches  zero.  Conversely,  very  near  the  origin,  f  approaches 
one  and  g  approaches  zero.  These  statements  may  be  made  in  light  of  the 
fact  that,  at  large  distances  from  the  origin  of  the  real  coordinate, 
the  magnetic  field  and  the  plasma  should  each  be  relatively  undisturbed 
if  only  a  finite  time  has  elapsed. 

With  this  much  knowledge  in  mind  about  the  behavior  of  the  variables, 
an  attempt  may  be  made  to  at  least  partially  understand  the  equations . 
Taking  Eq.  (38) 


Integrating 


zf '  ^  zfg* 
2g2  2g3 

(44) 

)•  -  -(  In  g 

)•  + 

Zffi* 

2ft  g3 

(45) 

f zfg'  dz 

J  2f'g3 

(46) 

[-/'?•/ 

^zfg'  dz  I 

2fg^  J 

(47) 

I 


1 
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Integrating  again 


Both  f  and  g'  go  to  zero  at  large  z;  thus  it  is  difficult  to  discuss 
the  behavior  of  the  right-hand  side  of  the  exponent  at  this  limit. 
Nevertheless^  it  might  be  expected  to  go  to  zeiro^  since  the  remaining 
terms  give  an  expression  for  f  similar  to  that  in  the  case  of  the  in¬ 
compressible  fluid.  This  is  vhat  might  be  expected  at  large  z  where 
g  is  nearly  undisturbed. 

Expressions  for  the  term  gf*  may  be  derived  from  both  equations. 
Eq.  (38)  may  be  written 

(sf)'  (|)'  (U9) 


which  Integrates  to 

8^'  “  -  /l  ^50) 

Also  from  Eq.  (39) 


z  f 


Both  of  these  may  be  integrated  to  give  expressions  for  f.  Unfor¬ 
tunately,  none  of  the  above  yield  any  new  information  about  the 
functional  form  of  either  f  or  g. 

\ 

It  turned  out  that  none  of  these  methods  were  really  helpful  in 
uncoupling  the  equations  to  obtain  an  expression  for  either  f  or  g. 
They  either  led  to  congplicated  integral  forms  or  approximations  too 
crude  to  be  useful.  It  is  possible  that  methods  for  treating  non¬ 
linear  differentiaQ.  equations  that  at  least  lead  to  further  simplifi¬ 
cations  in  this  case  might  exist.  Such  methods  are  discussed  in  texts 
like  Ince^^^  and  Kaplan^^^;  however,  their  application  is  sufficiently 
difficult  that  no  attempt  was  made  to  use  them  on  this  problem. 
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Another  possible  method  of  investigation  is  to  substitute  various 
well-knovn  functions  for  the  variables  and  see  if  they  bedance  the 
equations  at  asymptotic  values.  A  great  deal  of  time  and  effort  was 
spent  trying  this, , and  it  proved  vmequlvocally  fruitless .  The  actual 
functions,  thovigh  they  might,  be  well-behaved,  are  of  sufficient  com¬ 
plexity  that  no  closed  function  resembles  either  of  them.  This  is 
illustrated  by  noting  that,  in  the  greatly  restrictive  case  of  an 
incompressible  plasma,  f  still  had  the  form  of  an  open  Integral. 


V  Numerical  Integration 

In  order  to  obtain  specific  data  on  the  plasma-field  interaction, 
it  was  necessary  to  resort  to  numeriCEd  integration.  The  procedure 
used  is  described  in  Kaplan^^^  and  simply  consists  of  selecting  values 
of  the  functions  and  their  derivatives  at  a  point,  calculating  their 
values  at  a  point  close  by,  and  iterating.  Thus  the  method  is  a  numeri¬ 
cal  first  integration;  and  since  f  appears  to  the  second  order  in  the 
equations,  it  is  necessary  to  integrate  twice  to  obtain  values  for  it. 


The  equations 

are  written  in  the  form: 

f" 

.  .  - 

zf *  zfg* 

(54) 

8 

2g2  2g3 

f 

-  (f)' 

(55) 

a* 

_  2 

(56) 

0 

r 

z(l-f^)  ^ 

(i-f^) 

where  the  first  equation  is  Eq;  (38),  and  the  third  is  Eq.  (39).  In 
doing  the  calculation,  values  for  f  and  g  and, their  derivatives  were 
selected  at  large  positive  z  and  the  iteration  carried  toward  the 
origin.  These  initial  values  were  adjusted  until  f  approached  one  and 
g  approached  zero  at  the  origin  and  both  held  their  asymptotic  values 
at  large  z.  It  would  seem  more  accurate  to  start  at  the  origin  and  look 
for  asynptotic  values  of  the  functions  by  iterating  away  from  the  origin, 
but  it  was  extremely  difficult  to  estimate  the  first  derivatives  at  this 
point. 
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The  results  of  the  numerical  integration  are  plotted  in  Figure  3. 
The  variables  are  left  in  the  dimensionless  form  of  f  and  g  and  are 
plotted  as  functions  of  the  transformed  coordinates. 

These  transformed  coordinates  were  then  recast  in  the  form  of  the 


reail  coordinates  by  means  of  the  inverse  transformations. 


Ii 


dX 


(57) 


This  was  also  done  numerically.  The  variables  f  and  g  are  thus  also 
plotted  as  functions  of  x  and  t  in  Figure  4.  Because  the  boundary 
condition  x  — *  >  oo  as  X  — ►  0  is  awkward  for  numerical  ceilculation, 
the  origin  of  the  reel  axis  was  determined  as  follows: 

g  dx  (58) 


so  that  there  is  an  Xq  such  that 


It  was  then  assumed  that  Xq  and  x  «  0  correspond  to  g  ■  ^  because  the 
plasma  does  not  move  as  a  mass,  due  to  the  balance  between  the  magnetic 
and  gas  pressures.  This  yielded  Xq  ■  .35  from  Figure  3  Xq  =  .348 
from  a  graphical  integration  of  g  over  x  in  Figure  4.  These  figures 
thus  determine  the  magnetic  field  strength  and  the  plasma  density  in 
;x)th  space  and  time  as  the  diffusion  process  takes  place.  By  also 
calculating  temperature  and  entropy  changes  in  the  gas,  its  total  be¬ 
havior  may  be  considered  adequately  defined. 

The  temperature  of  the  gas  at  each  point  was  calculated  from  the 
x-axis  curves  of  f  and  g  as  follows:  since 


p  -  ^RT 


(60) 


then 


_E_ 


B 


B 


Changing  to  dimensionless  variables 

V  (1-f^) 


T 


8 


(61) 


a 


* 


i 


(62) 
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which  Is  plotted  In  Figure  5  with 


(63) 


Since  the  process  was  assumed  adiabatic,  the  specific  entropy  was 
calculated  from  the  x-axis  curves  of  g  and  f  by  using  the  expression 


8  -  80  +  Cp  In  +  °v  10  ^ 

e. 

Changing  to  dimensionless  variables 


2  2 
B  -  B 
o 


s  -  Sq  +  Cp  ln(|)  +  c^  In 
=  So  +  Cp  ln(i)  +  Cy  In  (1  -  f^) 


(64) 


(65) 


Or,  subtracting  Sq  and  dividing  by  Cy 

-  ln(i)  +  In  (1  -  f^)  (66) 

'-y  O 

Tlien  the  increase  in  entropy  as  a  function  of  x  and  t  is  given  by 

^  s  =  fimS  which  was  calculated  and  plotted  in  Figure  6  using 
y"  =  5/3.  The  total  Increase  in  entropy  as  a  function  of  time  was 
obtained  by  graphically  integrating  the  curve  in  Figure  6,  and  was 
found  to  be 


All  of  the  calculations  for  these  graphs  are  only  to  slide-rule 
accuracy.  Because  of  this,  and  because  of  thie  error  involved  in 
estimating  the  boundary  conditions  for  the  f  and  g  curves  in  Figure 
3,  the  curve  shapes  and  the  derived  quantities  may  only  be  considered 
to  be  correct  to  someidxat  better  than  an  order  of  magnitude.  Greater 
accuracy  would  require  the  services  of  a  conqputer. 
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